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Abstract 

OO , We consider the extension of the thermodynamic Bethe Ansatz (TBA) to cases in which 
additional terms involving higher conserved charges are added to the Hamiltonian, or in 

C/3 ' which a distinction is made between the Hamiltonian used for time evolution and that used 

Cu , for defining the density matrix. Writing down equations describing the saddle-point (pseudo- 

^*J' equilibrium) state of the infinite system, we prove the existence and uniqueness of solutions 

■4-J I provided simple requirements are met. We show how a knowledge of the saddle-point rapidity 

C . distribution is equivalent to that of all generalized inverse temperatures, and how the standard 

^ ' equilibrium equations for e.g. excitations are simply generalized. 

^ ; 1 Introduction 

Exactly solvable models of quantum mechanics are such that the basic Hamiltonian Hq enforcing 
'■^ ' time evolution of the theory is part of an extended family of conserved charges Q„ in involution, 

j^ . Hq, Qn = 0, 



^n7 Wr. 



0, V n, m. (1) 



Physically interesting cases correspond to charges expressed as integrals of (quasi-)local operators 
and the existence of this set ultimately enforces the notion of factorized scattering associated with 
integrability. Simple charges usually include a number operator N = Qo, the total momentum 
f^ ' operator P = Qi and obviously the Hamiltonian itself, Hq = Q2. The rest, Qn, n > 2 are not 

CO , usually explicitly invoked or used (although they are in principle available from the Algebraic 

Bethe Ansatz [T] via trace identities), partly due to the fact that expressions for higher conserved 
fT^ i charges quickly become unwieldy as one moves up their hierarchy. 

^D ' It is a completely trivial fact that any functional combination of conserved charges is itself a 

^N I conserved charge. Viewing the charges as the 'basis vectors' of the space of conserved charges, one 

can thus consider the generic expression 

"O' 7J({/3}) = ^/3„Qn (2) 

- - - as defining a generalized Hamiltonian depending on a (potentially infinite) set of real parameters 

/3„. The interest of such constructions is twofold. First, and most obviously, it permits to easily 
extend the 'phase space' of integrable Hamiltonians considered, by using the charges as represen- 
tatives of competing interactions [2j [3j HI [3 El |7] , most typically by including only one or a few 
additional charges as compared to the basic Hamiltonian Hq. Second, there exist situations in 
which one would like to consider functional averages (expectation values) using a density matrix 
which is not given by the usual Gibbs weight involving the Schrodinger Hamiltonian Hs used for 
time evolution, but rather with a more general density matrix defined as [5J E] 



Pgge 



^e-5:.^"Q" (3) 



which is now generally known as the generalized Gibbs ensemble (GGE), and which purports to 
represent the large-time density matrix resulting from a quench into an integrable system. In this 
case the generalized (inverse) temperatures /3„ are set by somehow solving the self-consistency 
problem for the initial conditions 

{Qm)^Tr{Q,ne'^-^"'^"}/ZGGE m = 0,l,2,... (4) 

in which Zgge — Tre^^"'^"'^" is the generalized partition function. Precisely how this is done 
depends on the situation; for free models it can be done explicitly by exploiting the factorization 
into separate momentum sectors (see for example [21 UHl [HI [HI [I3])- In the generic interacting 
case the trace in ^ cannot be taken explicitly. However, in the thermodynamic limit it is justified 
to consider the saddle-point of the partition function, which we will derive in the next section for 
the Lieb-Liniger model. On the other hand, defining the charges themselves and performing the 
averages in Q typically represent an insurmountable problem. 

It is our purpose here to provide a rather simple extension of the well-known thermodynamic 
Bethe Ansatz (TBA) to these situations, and to discuss the existence of solutions and the actual 
implementation of calculations in these cases. 

2 Generalized TBA for the Lieb-Liniger model 

For concreteness, we will illustrate the ideas by considering the Lieb-Liniger model, whose Hamil- 
tonian is 

H^= f dx{d^-^\x)d^-^{x)+c-^\x)-^\x)-^{x)-^{x)] . (5) 

Jo 

We restrict ourselves to the case of positive coupling constant c > 0, in which energies per unit 

length remain finite in the thermodynamic limit. As is well known, this model can be solved via 

the Bethe Ansatz. For a finite number of particles N in a periodic system of length L, eigenstates 

are fully characterized by a set of real rapidities Xj satisfying the Bethe equations 



zAv L 



Xj — Xi + ic 



i^, A. - ^^ - *c 

The eigenvalues of N, P and Hq on an eigenstate |{Aj}) are respectively Qq ^ N = X^i -^?' 

Qi — J2j ^j ^nd Q2 = J2j^^- The logic extends to all higher conserved charges Qn, whose 
eigenvalues are simply given by the power sum symmetric polynomials 

Qn\{X})^Qn\{X}), Q„ = ^AJ. (7) 

3 

The generalized Hamiltonian is therefore diagonalized according to 

00 Af JV 

H{{mw)^mmw), E{{i3}\{x})^Y.Y.f^^^^i^T.^o{Xj) (s) 

n=0 j=l j=l 

in which we have defined the function 



eo(A) = ^/3„A" (9) 



n=0 

by interpreting the coefficients /3„ as those of its power series. 

For a generic eigenstate in the thermodynamic limit, one can approximate the distribution of 
rapidities by a continuous density function p(A) taking positive (or zero) values for all A, the Bethe 
equations © transforming into an integral equation [14]: 

P(A)+P/.(A) = — +a2*p(A), a2iX) = - " (10) 

2tt tt A^ -|- c^ 



in which we use the standard convolution notation /*g(A) = /_ dX' f{X — \')g{\'). The function 
Ph{X) then represents the density of available quantum numbers (holes). For a given p(X) one can 
then easily compute the expectation values of the conserved charges as [I] 



dAA"p(A). (11) 

-OO 



poo 

L 

' — OO 



Following the arguments given by Yang and Yang [TS] , we now consider a generalized partition 
function 

Z^ I V[p\e-^^P'P'^^P^^ (12) 



(note that we have explicitly written ph as a functional of p by using the Bethe equations ©) in 
which the measure is given by a generalized Gibbs 'free energy' (although this is strictly speaking 
not an energy anymore, but a dimensionless quantity) functional 

G[p, Ph] = ^ PnQn - S[p, ph] (13) 

where the entropy is given to leading order in system size as [151 116) 

/CXD 
d\ [{p + Ph) \n{p + Ph) - plnp- Ph In ph] ■ (14) 

-OO 

In the thermodynamic limit we can evaluate the partition function in the saddle-point approxi- 
mation. This leads to the saddle-point condition 



PhjX) 
P(A) 

By proceeding as usual and defining the function 



In ^ = ^ /3„ A" - a2 * ln[l + p{X)/ph (A)] . (15) 



the saddle-point condition can be rewritten as 

£(A) + a2*ln(l + e-^(^)) = £o(A) (17) 

which we call the generalized TBA equation. This is a straightforward generalization of the usual 
TBA equations, the only difference being that the driving function eo{X) is given by the generic 
polynomial ^. One could even think of lifting the restriction of Eq to polynomial functions, and 
consider functions with isolated singularities. In view of the applications we have in mind, we will 
however not consider these more general cases here. 

The equilibrium (i.e. saddle-point) state of the generalized distribution (I13p is completely 
determined by ((TU)) and P7)) . One can rewrite (ITUl) to eliminate ph{X) obtaining 

p(A)=^?(A)(^i-+a2*p(A)) (18) 

in which ^(A) is as usual called the Fermi weight and is defined as 

1?(A) = ,^/^^^ ,^, = ^—-r. (19) 

^ ■' p(A) + p,,(A) 1 + eM^) ^ ' 

As in the usual case, many expectation values only depend on p(A) and ^(A) (see for instance 
[T71 [T51 [HJ [201 m]); these results can therefore automatically be generahzed to solutions of the 
generalized TBA equations. 



For a given set of /?„ the saddle-point state is completely determined by equations (fTT)) and 
T8l) . Conversely, for a given e(A) the Lagrange multipliers can be determined explicitly from ((T7)) . 



/3«= ^(e(A) + a2*ln[l + e--W]) 



(20) 



A=0 



One can proceed in either of two ways, depending on the available data. Either for a given 
set f3n one solves for e(A) using (fT7)) and from there finds p(A) using (fT5|) . This would be the 
procedure to follow for example in the case where the explicit values of the f3n correspond to 
specific 'user-defined' perturbations to the original Hamiltonian [2l [3l [7] , or to situations in which 
(by some miracle) the generalized inverse temperatures of the GGE ensemble are known. 

A dual interpretation and use of the generalized TBA equations consists, in fact, in starting 
from a given p{X), and solving for e(A). The Lagrange multipliers, which ultimately encode all 
characteristics of the saddle-point state, can then be computed via ([^0)1 . While it might seem 
strange to expect the distribution p{X) to be given as input, this case is in fact the one which 
occurs when one uses (numerical) renormalization around a Bethe Ansatz-solvable point [22j . This 
case in fact represents possibly the most immediately useful application of the gTBA equations, 
since the knowledge of p(A) then allows to explicitly compute GGE predictions without explicit 
knowledge of the set {Pj}, by using the generahzed TBA equations to relate p(A) to the physical 
distribution e(A) actually used in the averaging. In the case of the Lieb-Liniger model, in fact, 
the problems with attempting a direct implementation of the GGE according to the prescription 
discussed in the introduction, are compounded by the fact that the conserved charges are not 
properly normal-ordered objects [231 IM] whose expectation values can easily be computed using 
standard methods. The alternate route which is offered by the generalized TBA equations is thus 
actually the only practical one available. We refer the reader to [22 for a parallel exposition of 
the gTBA approach together with its explicit implementation in specific quench cases. 



3 Solving the generalized TBA equations 

We can show that a solution of the generalized TBA equation ([T7)) exists and can be found via 
iteration, in complete parallel to the traditional case, provided that two very simple conditions 
are fulfilled. The bare energy 

eo(A)=^/3„A" (21) 

n 

should be a) bounded from below, and b) be such that limA^±oo ^oW = -|-cxd. 

The proof we give is essentially a repetition of the traditional one offered in [IS]. Let us 
construct the following sequence of functions 



where 



en+i(A) = eo{X) + A[en{X)] 



A[e{\)] = - I dX'a2{\ - A') In (l + e'^^^'^) 



(22) 



(23) 



The proof contains two steps. First we will show that for every A, the sequence of functions is 
strictly decreasing, 

£o(A) > ei(A) > . . . > e„(A) > £„+i(A) > . . . (24) 

Secondly, we will show that this sequence is bounded from below so that the limit 



e[X) = lim £„(A) 



(25) 



exists and is a solution of (fT71) . 

From ([231) we see that A[e(A)] < for all e(A). In order to show that yl[£„+i(A)] < yl[e„(A)] 
we consider 

M[£„(A)]=/ dA'a2(A-A')- -jrj.5en{X'). (26) 



For fe„(A) < we have (5^[e„(A)] < hence we arrive at ([M| . To prove that e(A) is bounded 
from below is more comphcated. First we speciaUze to cases in which So{X) is symmetric in A and 
is monotonically increasing for positive A. By considering 

+ dXa2{X~X)— —-—.—-— - (27) 



dX dX J_^ ^' 'l + e^"(^') dX' 

we can prove by induction that en(A) is also symmetric in A and monotonically increasing for 
positive A. From this follows that A[e„(A)] is monotonically increasing as a function of A, and 
goes to zero in the limit A — )■ oo. Using this fact we can write the following inequality 

£„(A) = eo(A) + ^[e„-i(A)] > eo(A) - eo(0) + e„(0) (28) 

where the equality holds for A = 0. Using this inequality we can write another inequality using 
(I22ll and (1261) 



en+i(0) > £o(0) - / dX'a2{0 - A') ln(l + e-(^°(^')-^°(o)+""(o»). (29) 

Next, we define the function 

/oo 
dX'a2{0 - A') ln(l + e-(so(V)-eo(o)+x)) (gg) 

-oc 

roc 

^eo{0) + x- / dA'a2(0-A')ln(e" + e-(^°(^')-"°(°») (31) 



so that (P51) can be written as £„+i(0) > /(e„(0)). The function f{x) increases monotonically and 
is bounded from above, f{x) < eo(0). One can also show that f{x) — x decreases monotonically 
and that its image is (—00,00). Thus the equation f{xo) — Xq = has a unique solution. For a 
given eo(A) we can now determine xq from the equation 

eo(0) = / dX'a2{0 - A') ln(e^° + e-^'''^^'^^'^^^^). (32) 



We shall now prove by induction that e„(0) > xq for all n. First one notes that eo(0) > f{xo) — xq- 
Next, suppose that e„(0) > xq. One then has from e„+i(0) > /(e„(0)) and the monotonicity of 
f{x) that 

£n+i(0)>/(£„(0))>/(xo)=a;o. (33) 

This proves the inequality 

£«(A)>eo(A)-£o(0) + a:o V A. (34) 

Combining this with (|24)) we have proved that the solution of (J17p can be found by iteration for 
a symmetric £o(A) which increases monotonically for positive A. Consider now a general £o(A) 
not satisfying these conditions. Since there obviously exists an £o(A) that is symmetric in A, 
increases monotonically for positive A and £o(A) < £o(A) for all A, then from (pS)) follows that 
A[£„(A)] < A[£„(A)] for all n and A. Hence, £(A) is bounded from below as well. 

4 Uniqueness of the solution 

We can also show that a solution of (IT71) extremizes the generalized free energy, following the lines 
of [TS]. We will also show that this solution is unique. 

Consider two solutions of the Bethe equations pi and p2- It is clear that xpi + {1 — x)p2 for 
< X < 1 is also a solution. Using this property we can define an action X{L, {/3„}, p) by 



^00 
X = L I dX 

-00 



p{X) ^ /3„A" + p\np + ph \nph - {p + Ph) ln(p + Ph) 



(35) 



and vary it with respect to p. Consider p — po + xpi where po and pi are two independent sohitions 
satisfying (fTO|) . The variable x takes real values. We can differentiate the action X with respect 
to X 



-— = L / dXpi ^ /3„A" - e(A) - / dA'a2(A - A') ln(l + e^^^^')) 



(36) 



Next, using (IT^ . we can compute 

de _l + e~^ f 1 
dx p 

Now, 



(y^/"^^'«2(A-A')pi(A')-Pi). (37) 



By first performing the integral over A and then using (j37p this can be simplified as 



dx'^ J_^ \ dx J 1 + e-^(^) 

Since this is true for any p we conclude that X is convex, hence it has a unique minimum. 

5 Thermodynamics 

Now that we have established saddle-point (equilibrium) state we can write down the standard 
thermodynamic relations. First we write generalized free energy (J13p as 



G^L p(A)(£o(A) - £(A)) - (p(A) + ph{X)) ln(l + e-'^^^)d\. (40) 

J —oo 

At the saddle-point state we use P7)) and ([TS]) to further simplify the expression 

G = L \n{l + e-''^^'>)dX. (41) 

27r y_oo 

We can easily show that G satisfies the following thermodynamic identities corresponding to the 

Hellmann-Feynman theorem 

f)G 

We first differentiate ([TT]) with respect to /?„ 

^'^^^=X-+ra,iX~pnp)^dp. (43) 



d(3n J-oo dp, 

Multiplying both sides with p[X) and integrating over A gives 

^p(A)dA ^ £ A>(A).A + £ {pip)/m - i^) ^(.)^c^.. (44) 

We have used (fT5| in the second term on the right hand side to eliminate the integral over A. 
From this follows the identity 

dG L f°° 1 9e(A)^^ ,- , ,^^, 

- ' ' 'dA= (g„). (45) 



9/3„ 27r 7_^ 1 + e^(^) dPn 
For the expectation values (Qn) we can also define corresponding generalized susceptibilities as 

d^G 



dPl 



{{QD - {Qn?) ■ (46) 



Using similar tricks as before we obtain 

which is essentially identical to ([5^ . In summary, we can simply state that all the usual equilibrium 
thermodynamic equations remain unchanged as compared to the usual case, provided one uses the 
generalized £(A) and p(A) functions. 

6 Excitations 

In order to study excitations upon the saddle-point state it is instructive to first go back to a state 
with a finite number N of particles. To a state with quantum numbers Ij corresponds a set of 
rapidities Xj that satisfy the Bethe equations 

N 

A,L = 2^/,-^0(A,-A,). (48) 

1=1 

Now we consider an 'excited' state with quantum numbers /' and the rapidities A'- satisfying 
different Bethe equations 

N 

XrL^2n^-Y,0{\',^Xd (49) 

i=l 

such that Ij = /,' except when j — n. We make the usual assumption that for all j ^ n the 
difference between Xj and A' is small. One can introduce a shift function x(-^) 

(A;-A,)L = x(A,) forjV«- (50) 

Going back to the thermodynamic limit one can show using the arguments of Lieb [33] and Yang 
and Yang ^15] that x(A) is determined by the following integral equation 

/oo 
a^ix - ^)(x(A) - x{^A)p{^J)d^l + e[x - a„) - ^(a - a;). (51) 

-CXj 

Writing the back-flow as g[X) ~ x(A)(p(A) + Ph[X)) we obtain 

g{X) = j a2{X-ii)gif,Wp)dti + — {0iX~X,,)-eiX-X'J). (52) 

The momentum difference and energy (as measured from the basic Lieb-Liniger Hamiltonian) 
difference between the two states are 

/oo 
g{X){){X)dX, (53) 

/OO 
2Xg{X){){X)dX. (54) 

-oo 

In the case in which the density matrix Hamiltonian does not coincide with the Schrodinger 
Hamiltonian, we have that AE ^ e(A^) — £(A„), therefore e(A) cannot be interpreted as the energy 
of fundamental excitations anymore, as opposed to the usual case. However, if the Schrodinger 
Hamiltonian is by definition the generalized Hamiltonian H — Heff — X^n/^nQm i^^ excitations 
are 

AE^jf = 5](eo(A;.) - eo(A,)) = soiK,) - eo(A„) + J^^ ^^g(X)HX)dX (55) 

= e(A;)-£(A„). (56) 



Hence e(A) are the excitations of the effective Hamiltonian. As in the usual case, it is straightfor- 
ward to prove that a finite number of simultaneous excitations is simply the sum of the individual 
elementary excitations 

APiiXp^ }, {A,^ })^J2 ^P(^P. ' ^'^. ) (57) 

i 

AS({ApJ,{A;.J) = ^Ai^(Ap^.,A;.,), (58) 

j 

this relation remaining true provided the density of excitations thus created remains zero in the 
thermodynamic limit. The whole bulk of knowledge about correlations of the equilibrium Lieb- 
Liniger gas JLj can thus be easily adapted to the generalized cases. 

7 Conclusion 

Although we have concentrated on the Lieb-Liniger model, similar reasonings are (almost trivially) 
transportable to other models. The generalization of generic TBA equations (with strings, nested, 
etc.) is straightforward since only the basic driving function (bare energy) needs to be modified. 
In cases other than the repulsive Lieb-Liniger model however, proving existence and uniqueness 
becomes even more intractable than in the traditional case, although explicit solutions can easily 
be found in practice. 

In summary, we have provided a simple extension of the traditional TBA equations to cases in 
which higher conserved charges come into play, either via extended Hamiltonians or via density 
matrix averaging differing from the usual Gibbs procedure. For the Lieb-Liniger, existence and 
uniqueness is proven provided the driving function fulfills very simple conditions. The use of the 
generalized TBA equations provides a practical road to the actual implementation of generalized 
Gibbs ensemble averages, or other similar ones, by offering an alternate route which obviates the 
need to work with the actual expressions for the conserved charges and their associated inverse 
temperatures. 
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